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The scattering of conduction electrons by crystalline electric field (CEF) excitations may enhance 
their effective quasiparticle mass similar to scattering from phonons. A wellknown example is Pr 
metal where the isotropic exchange scattering from inelastic singlet-singlet excitations causes the 
mass enhancement. An analogous mechanism may be at work in the skutterudite compounds 
Pri_2,La3;Os4Sbi2 where close to x=l the compound develops heavy quasiparticles with a large 
specific heat 7. There the low lying CEF states are singlet ground state and a triplet at A = 8 
K. Due to the tetrahedral CEF the main scattering mechanism must be the aspherical Coulomb 
scattering. We derive the expression for mass enhancement in this model including also the case 
of dispersive excitations. We show that for small to moderate dispersion there is a strongly field 
dependent mass enhancement due to the field induced triplet splitting. It is suggested that this 
effect may be seen in Pri_a;LaiOs4Sbi2 with suitably large x when the dispersion is small. 

PACS numbers: 71.27.+a, 75.10.Dg, 71.38.Cn 



I. INTRODUCTION 



The fiUed-skutterudite compound PrOs4Sbi2 has recently obtained considerable attention. There are several reasons 
for that. It is a heavy fermion (7 ~ 350 — bOOmJ / molK^) superconductor with a transition temperature of T'c(Pr) = 
1.85 K. This temperature is larger than the one of the related system LaOs4Sbi2 which is Tc(La) = 0.74 K. A number 
of experiments, like those on Sb-NMR relaxation rate in Pri_3:Laa;Os4Sbi2 pLj] suggest that the superconducting 
order parameter is of the conventional isotropic s-wave type with possible admixture of higher harmonics depending 
on the Pr content [2*1 . However, questions and ambiguities remain. They concern the experimental findings for the 
penetration depth [3,, a-nd initial studies of the thermal conductivity in a rotating magnetic field ^5]. For example, 
the former suggest a possible nodal structure while the latter in addition points towards two distinct superconducting 
phases. The observed enhancement of the superconducting transition temperature by more than a factor of two when 
La is replaced by Pr seems surprising at first sight. It is well known that when Pr ions are added as impurities to an s 
wave superconductor like LaPbs it suppresses the superconducting transition temperature rather efficiently. So why 
does it enhance Tc in the present case? Since the phonons in LaOs4Sbi2 and PrOs4Sbi2 are very nearly the same, the 
enhancement must come from the two 4/ electrons which Pr'^^ has. The heavy fermion behavior of PrOs4Sbi2 seems 
puzzling too. It shows up in a large specific heat jump AC/Tc ~ 500 mJ/(mol K^) at Tc and also in a large effective 
mass in de Haas-van Alphen experiments. The Kondo effect cannot be the origin of the heavy quasiparticles since the 
4/^ electrons are well localized with a Hund's rule total angular momentum J — A and a non-Kramers ground state. 
The key to the enhancement of Tc and the formation of heavy quasiparticle excitations lies in the crystalline electric 
field (CEF) splitting of the J = 4 multiplet (Sect. EJ, together with the aspherical Coulomb scattering mechanism 
of conduction electrons from CEF excitations [1, 01^ The self energy and effective mass enhancement due to this 
mechanism will be calculated in Sect. Hill The CEF states, their excitation energies and matrix elements are modified 
by an external field. The ensuing effective mass dependence on the field which is the main topic of the present work 
is calculated in Sect. llVl for dispersionless excitations and in Sect. |V] for the case with dispersive quadrupolar excitons 
. Some numerical results are discussed in Sect. IVII and Sect. IVIII finally gives the conclusions. 



II. THE CEF STATES OF Pr IN T^ SYMMETRY AND THEIR INTERACTIONS 



From inelastic neutron scattering the CEF energy levels are known. The compound has tetrahedral Th site symmetry 
for Pr. The data are explained best by a CEF ground state 



FIG. 1: Left (a): Fermi surface of n.n.n. tight binding model in hole representation in the bcc Brillouin zone. Right (b): 
schematic Fermi surface in electron representation in the 2D projected Brillouin zone. It consists of spheroids around the 
equivalent H-points 0, 0). The polar angle ^ of q is given hy d = ^{ir — 9'). Furthermore we have q — 2pF sin ^ = 2pF cos 6 
where pp is the Fermi momentum. The geometric restrictions require < g < 2pF and < ^ < |-. 

with a low- lying triplet excited state at an energy of A = 8K d, [lol| . The other CEF levels are so high in energy 
that they can be neglected. The Ft triplet state of Th symmetry is a superposition of two triplets r4 and Ts of Oh 
symmetry. More specifically one finds [1, [lo| 



I Tt,m) = ^l-d^\ r5,m) + d| r4,m) , m = T..3 (2) 
with states of Oh symmetry given by 

|r5,±) = ±y||±3)TY^|Tl) |r5,0) = yi(|+2)-|-2)) 

\T,,±) ^ T^\t3)T^\±1) |r4,0) = yi(|+4)-|-4)) . (3) 

The conduction electrons interact with the CEF energy levels of the Pr^+ ions. The most important ones are the 
isotropic exchange interactions and the aspherical Coulomb scattering. The former is of the form 



Hcx{i) = -2 {gj - 1) Jcx ^ (s<T<T' • Ji)4-q<T''^k<T (4) 

where c|j.^(ck(j) are the creation (annihilation) operators for conduction electron with momentum k and spin a while 
s is their spin operator. Furthermore gj is the Lande factor. The aspherical Coulomb interaction in local orbital basis 
is given by pTj 
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FIG. 2; Static (normalised) quadrupolar susceptibility Axq as function of magnetic field with (d^ — 0.067 or |d| — 0.26) and 
without (d^ — 0) tetrahedral CEF. Full lines correspond to evaluation with Eq. (|24|l while the dashed line corresponds to 
the low field expansion Eq. (|26|l . Approaching the Fa — F^ level crossing leads to an increasing XQ Sund mass enhancement 
which becomes singular at the crossing he/ A — (1 — 5'^)~^ . The latter is pushed to higher field for increasing which reduces 
the increase of XQ- For > 0.42 XQ decreases with field strength because the tetrahedral CEF leads to Fs — F? repulsion. 
The mass enhancement is proportional to the quadrupolar susceptibility with 5m* /m — ^^//(Axq) and g^ff — {g/A)f. For 
d^ = 0.067 and g^ff = 0.077 (Sect.E]) one has {Sm* /m)h^o = 16. 

Here Ckima- destroys a conduction electron with momentum A: — |k|, in a / = 2 state with azimuthal quantum number 
m and spin a and c\.,^^ creates one with momentum fc' in a Z = state. The Coulomb integral I2 is defined, e.g., in 
Ref. [Til and Q2 is the quadrupolc moment of the Pr'^+ ion. The operators Y{^{3) are given by 

y^" = (2/3)^/'[3j2- J(J+l)]/iVj 
K,* = ±{J,J^ + J^J,) /Nj 

- {J^f/Nj (6) 

with Nj = (2/3) "'^/^ (2 — J). Hac leads to a transfer of angular momentum I = 2 between the conduction electrons 
and the incomplete 4/ shell. It is a quadrupolar type of interaction. 

An important feature of PrOs4Sbi2 is the experimental finding that the low-energy triplet state has a small value of 
\d\ = 0.26 with the implication that the inelastic scattering of the conduction electrons is predominantly of quadrupolar 
character. With this information the two features pointed out above, i.e., the increase of Tc when La is replaced by Pr 
and the heavy quasiparticle mass can be understood quantitatively Q . As has been known for a long time quadrupolar 
inelastic scattering of conduction electrons by low-energy CEF levels enhances Cooper pairing since these excitations 
act similarly as a localized phonon mode. The difference is that phonons are related to changes in the ion position 
while intra-atomic quadrupolar CEF excitations are related to changes of the 4/^ wavefunction. 

Also the heavy quasiparticle mass is related to the inelastic scattering processes of the conduction electrons. This 
feature has been previousl y ex ploited to explain the mass enhancement found in Pr metal using the isotropic dipolar 
exchange interaction Hex UM- mentioned above in Pri_a;La2;Os4Sbi2 the aspherical Coulomb scattering Hac is 
dominant over exchange scattering. For this model a quantitative calculation of the changes in Tc and the mass 
enhancement as function of Pr concentration are found in Ref. 0- 

The aim of the present communication is to extend the previous calculations by including the effect of an external 
magnetic field on the mass enhancement. A field splits the triplet states and leads to a decrease of the excitation 
energy of one of the three states, at least for small tetrahedral admixture d as in PrOs4Sbi2. Therefore an increase 
of the effective mass with increasing field is expected in this case. 

The situation is different for Pr metal mentioned before where ground state and first excited state are two singlets. 
There a magnetic field pushes the two energy levels apart and hence increases the excitation energy. As a consequence 
the effective quasiparticle mass decreases with increasing external field in agreement with experimental findings . 
In the present singlet-triplet model this case would be reahsed for \d\ > 0.65. 
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III. SELF ENERGY AND MASS RENORMALISATION 

We start out by specifying the electronic part of the Hamihonian for the system Lai_a:Pra;Os4Sbi2. It is of the type 

H = Hc\ + i?cEF + Hp^c + ^^cx + Hz ■ (7) 

Here Hd is of the conventional form and need not be explicitely written down. It contains the conduction band 
dispersion which may be described by a n.n.n. tight binding model (Tsj according to 

^ko- = t cos -/cj; cos -fcy cos i^^z ^ t'(cos fcj. + cos ky + cos kz) (8) 

with t=174 meV and t'=-27.84 meV. The transfer integrals t and t' are chosen so as to reproduce the observed 
linear specific heat coefhcient 7 — 36mJ/(mole K^) of the non-f reference compound LaOsi2Sbi2 p^. In the electron 
picture the associated Fermi surface consists of H-centered spheroids with a Fermi radius pp ~ 0.7(27r/a), see Fig.[Tl3. 
Aside from subtle effects this is quite similar to the LDA FS in Ref. 17 (note that in this reference the FS is depicted 
in the hole picture as in Fig. [T^). It corresponds to a single band originating in Sb-4p states. 

The CEF and Zeeman Hamiltonian are 

i?CEF + ifz = ^^r|r„(z))(r„(i)|+gjAiB^J(*)-H . (9) 

The external magnetic field is denoted by H, gj is the Lande factor and fiB is Bohr's magneton. Furthermore i labels 
the Pr sites, and |r„) denotes the singlet ground-state jr^) and the triplet jF") (see Eq. ^) with energies Eg = 
and £"4 = 8 K, respectively. We assume that not only the phonons but also their local interactions with the electrons 
are independent of partial replacements of La by Pr. 

As pointed out before the system has Th symmetry but since the CEF transition can be reduced to those between 
states of cubic symmetry (see piBI) ) we specialize ([5]) to cubic symmetry. In that case the aspherical Coulomb 
interaction written in a basis of Bloch states becomes 

H^c(^)^9Y. E OL;39a9>ctq,Ck.e^'^^' (10) 

with Oap = -^{JaJp + JpJa), Qa — la/l'il and a/3 — yz,zx,xy denoting the three quadrupole operators with Fs 
symmetry. The remaining F3 quadrupole terms are neglected since they do not couple to the singlet-triplet excitations. 
The coupling constant g may in principle be determined by experiments. A derivation of (jlOp may be obtained from 
Ref. M 

In order to determine the effective mass m* of the quasiparticles at zero temperature one must calculate the Green's 
function of conduction electrons. 

G(p,w) = ^ . (11) 

uj - e(p) - S(p,w) 

The effective mass enhancement due to interactions of the conduction electrons follows from 



= 1 - 



duj 



(12) 



where pf is the Fermi momentum and m is the reference quasiparticle mass including band effects and effects of 
electron-phonon coupling. Neglecting vertex corrections the irreducible electron self energy I](p,w) due to Hac is 
given by 

E(p,c.)=5'E/7CTl^S/3(q)l' /^^„(q,c.)G'(p-q,c.-c^') . (13) 

a/3,n 
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Here I?„(k, denotes the boson propagator of CEF excitations. It is related to the dynamical quadrupolar suscep- 
tibility of the CEF system. In the present case we will neglect effective RKKY type interactions between CEF states 
on different sites therefore the boson propagator is local (q- independent). The momentum dependence of the bare 
vertex A^o(q) is due to the quadrupolar form factors in Eq. (fTU|) . It is defined as 



Kf,{ci)^9qaqp{rs\o^p\ri^) . (14) 

The self energy due to exchange scattering which involves the magnetic susceptibility can be safely neglected because 
of the smallness of (P (see ([2])) as a more detailed investigation including matrix elements and coupling constants 
shows. The propagator of the local singlet-triplet boson excitations is given by 



D„(q,c.) = A.(c.) = T5^ . (15) 

Here the field dependent singlet-triplet excitation energies are given by Sn{H) — e"(i/) — £s{H) (n = -|-,0 -) 

The self-energy in Eq. can be evaluated following Migdal's integration procedure (see. e. g. [iMl)- This 
method exploits the fact that the summation over fermionic states in the vicinity of the Fermi surface separates into 
independent summations over energy and degeneracy variables. Since the relevant excitation energies are of the order 
of the CEF excitation 6n (H) the dominant contribution to the integral in Eq. (fT3|) comes from electronic states with 
|e(p)| ~ (5„(0) <C W. First one replaces the integral over q by an equivalent integration over p' = p — q where the 
external momentum p is kept fixed. Then one restricts the frequency integration to a shell |cij'| <^ ec around the Fermi 
surface such that i5q <C Cc ^ is fulfilled. Here 2W is the conduction band width. In this shell we may approximate 
the momentum space integral by 
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where 9' , cj)' are the polar and azimuthal angles of p'. Furthermore N{0) is the conduction electron density of states 
per spin at the Fermi energy (ep = 0). Using the analytical properties of the self-energy S(a;) which imply that 
sign lml]{uj) = — signw the e'-integration gives for 5n{H) 

E(.) .,^A.(0) E(.^'.|>l(r.|0.,|rr)P/ g (,^,,)?:,.^,^ W2arctan^^ . (17) 

We solve the self-consistency equation for the imaginary part ImI](cL;) 

ImS(a,) = ~g^N{0)Y,{qlql)\{T,\Oo.p\ri)\'' ^ sign(u; + p5n{H)) arctan (18) 

ap.n p=±l ' ^ h' n\ J)\ 

from which we subsequently deduce the real part by Kramers-Kronig transformation. The explicit form immediately 
shows that 

|ImE(c.)|<ViV(0)5](g^g|)|(r,|O„^|rnP = J . (19) 

a/3,n 

It is important to note that ImE(cj) exhibits discontinuous jumps at the energies corresponding to the singlet-triplet 
excitations. This feature is a direct consequence of the assumption that the CEF excitations are long-lived and 
dispersionless bosonic excitations. Of particular interest is the discontinuity at 5^{H) 

|ImE((5+(i?) +77) -IniI]((5+(iJ) +?7)| > --arctan^ . (20) 

This discontinuity in the imaginary part inevitably implies a logarithmic singularity in the real part ReE(aj) which, 
in turn, leads to an unphysical divergence in the effective mass for 5_|_(_ff) — > 0. 

In the limit e,, ^ 00 where 2 arctan -t — — - — )■ vr the result agrees with that of non-sclfconsistent second order 

perturbation theory. In this case differentiating the self energy with respect to oj under the integral and using 
integration by parts one finally gets from Eq. (|12p : 
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TABLE I: Singlet-triplet CEF states, levels and excitation energies in a magnetic field H || [001]. Here S„{H) = E"{H)~Eg{H). 
The eigenstates are given in terms of zero-field singlet-triplet states |0, 0) and ll,±), |1,0) respectively {h — q^bH). 
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— = l + g^N{0)fxQ{H)- 
m 

The directional average (with respect to polar and azimuthal angles 9, of q) for quadrupolar form factors 
/ — {QaQp) — TS '^^ ^ constant. Furthermore xq{H) in Eq. pT|) is the field-dependent static uniform quadrupo- 
lar susceptibility. Note that the form factor average can be trivially factored out as a constant (1/15) only because 
in the present local approximation the boson propagator is momentum independent, i.e. the CEF excitations are 
dispersionless. The more general case will be discussed below. 



IV. FIELD DEPENDENCE OF THE EFFECTIVE MASS: DISPERSIONLESS MODEL 

When a magnetic field is applied to the sample the field dependence of the effective mass is completely determined 
by that of the quadrupolar susceptibility in Eq. (j2ip . To calculate this quantity one first has to know the singlet-triplet 
excitation energies 5n{H) and the eigenstates and matrix elements in applied field. They were given by Shiina et al 
[9l.[Toj in closed form for field applied along cubic symmetry directions. We use these results in the following. The CEF 
eigenstates are denoted as singlet jFs) and triplet IF") (n — -l-,0,-), respectively. The CEF and Zeeman Hamiltonian 
can be easily mapped to a pseudospin basis [9|, llO| and then diagonalised. In pscudospin basis the zero-field singlet 
is denoted by |0, 0) and the triplet by |l,m) (m=l,0,-l). For field H || [001] the field-split CEF levels and mixed 
eigenstates are then given in Table HI The field dependence of Sn{H) has recently been determined by INS eperiments 

The mixing coefficients u, v are determined by the matrix elements of the dipolar operator J in the Zeeman term 
which may be expressed by a = 5/2 — 2cP, (3 — 2^/5/3^, 5 = (3/a. They are given by 

u = -S5n(y)[i(l - ^)]5; m = (1 - w^)^ 

A = [A^ +4(52/i2]5; h = gfiBaH (22) 

Note that a finite mixing v occurs only due to the tetrahedral CEF contribution {y 0) which leads to o? 7^ in 
Eq. ([2]) and hence (5^0. When d = {6 = 0) there is no mixing between [0, 0) and [1, 0) and consequently the energies 
of \Ts) and [Fj) will be independent of the field H. For nonzero d and v these two levels will repel with increasing 
field H . The other two triplet states [F^) have a linear Zeeman splitting of 2h independent of d. For S > the singlet 
ground state level Eg and lowest triplet level cross at a critical field he = A/(l — 5^) meaning 6+{hc) — 0. 

For evaluation of the effective mass we need the quadrupolar matrix elements in Eq. (|2ip . They may all be expressed 



in terms of the irreducible zero field matrix elements a' = ^(13- 20^^) and (3' = ^35(1 - d^). With their help and 
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defining 6' = a' / (}' one obtains the following nonzero matrix elements 



|(F,[0,,[F+)[^ = [(F,[0,.[[Fr)|^ = lf5'\u-6'vr = \mQ\^ 

|(F,[0,,[Fr)|^ - |(F,[0,.[[F+)[2 = ^p'\u + 6'vr = \m+\^ (23) 

0\|2 _ a/2 _ 1^0 |2 



|(F,[0.,[F^)[2 = /?' 




0.4 0.6 
(q,0,0) (r.l.u.) 




0.4 0.6 
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FIG. 3; Left (a): Dispersion of lowest triplet quadrupolar exciton uj+{q) at the critical field hc{gQ) where tj+(Q) becomes 
soft (Q in units of 2n/a). As gg is reduced the dispersion becomes flat increasing the phase space for low energy conduction 
electron scattering. Right (b): Field dependence of AxQ ^ 5m* /m for various strengths of intersite quadrupolar coupling 
qq (left). For small gg mass renormalisation close to he is large due to flat w+(q) dispersion. For larger gg the dispersion 
becomes stronger and Lj+{q) softens only in the vicinity of Q = (0,0, 1) leading to a much smaller mass enhancement at he- 
The curves correspond to gg given in the legend in decreasing order. For the value gg = 0.3 corresponding to PrOs4Sbi2 little 
field dependence of 5m* /m remains. 

Inserting the matrix elements and excitation energies in Eq. (PT|) and using Irrigp + |mQp| — — — 5'^)] 
we finally obtain th.e expression 



xq{h) 



2/3' 
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A'2 - W 



(24) 



where we defined A' = i(A + A). Without tetragonal CEF {(P = 0) we have A' = A = A and then the above 
expression reduces to 



xq{h) 



2^ 
A 



2A2 
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[l-t;2(l-^'2)] + l 



(25) 



For small fields {h <C A) the general Xq{H) in Eq. (P5)l may be expanded with a leading term ~ (/i/A)^ according 



to 



(26) 



The zero field mass enhancement without tetragonal CEF { cP = Q) is then obtained from Eq. (|21l25p simply as 



^ = 1 + .927V(0)3/^ 
m A 



(27) 



The states and energies in Table U are nominally derived for H || [001]. However it was shown in Rcfs. [1, [13 
that for low fields {h < A) they are the same for all field directions, i.e. approximately isotropic. Therefore the 
quadrupolar susceptibility derived above and the related mass enhancement are also approximately isotropic as long 
as h is appreciably below A. This condition is required anyway in the dispersionless case where the calculation is only 
valid for moderate fields when S+{h) is still large enough. 



V. INFLUENCE OF QUADRUPOLE EXCITON DISPERSION ON THE MASS ENHANCEMENT 



In the previous section we investigated a model of noninteracting local singlet-triplet quadrupole excitations. For 
appreciably large concentration of Pr ions in the system Pri_a;La2;Os4Sbi2 this is no longer justified. Due to effective 
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FIG. 4: Left (a): Variation with magnetic field of 5m* /m calculated from the self consistent calculation for finite band width 
~ 205„(0) (full line) and ec ^ oo (broken line). Right (b): Variation with magnetic field of the specific heat coefficient 7 
calculated selfconsistently for ec = 20S„{0) 



interactions between the 4f states on different sites the singlet-triplet excitations at A acquire a dispersion. Formally 
this is already included in the self energy of Eq. (jl3|) provided the boson propagator for a dispersive mode is used by 
replacing Sn — > w„(q) in Eq. (|15p . The dispersion is due to quadrupolar RKKY-type intersite interactions which are 
obtained in second order perturbation theory from TIac and given by [12.] 



^^0 = E^Q(*J)OW-0(-?) (28) 

(ij) 

where O = {Oyz,OzxTOxy) is the Fs type quadrupole. The sum is restricted to nearest neighbors and Kq is the 
effective quadrupolar coupling constant. This interaction leads to the field induced antiferroquadrupolar order from 
which the F5 symmetry has been infered [§] . Formally Hg may be obtained from an RKKY type mechanism in order 
^ in the coupling constant of Hac- In practice the n.n. term Kq is determined from the experimentally observed 
dispersion of the quadrupolar excitons tj„ (q) which are fully degenerate for zero field [2l[ . In finite field the dispersive 
excitation branches arc obtained by replacing the Hamiltonian in Eq. ^ with Hqef + Hq + Hz- Using a generalised 
Holstein-Primakoff approximation [l^ the three quadrupolar exciton modes at moderate fields are described by 



w±(q) = w(q)=F/i; ^o{^)=^{^) (29) 
where the zero-field dispersion is given by 



c.(q) = - Bl 

= A + zi^Q7q; Sq = -zifQ7q (30) 
111 

7q = COS -qx COS -qy COS -qz (31) 

Here z=8 is the coordination number and 7q the structure function of the bcc cubic lattice of Pr ions with 
momentum q measured in r.l.u. (27r/a). The width of the exciton dispersion is controled by the effective quadrupolar 
coupling constant Kq or in dimensionless form by qq = z^'"^ Kq/ ^. From the analysis of the AFQ phase diagram 
[l^ and experimental zero- field dispersion [2]| one may deduce qq ~ 0.3 in PrOs4Sbi2. The minimum of the 
dispersion occurs at the bcc zone boundary wave vector Q — (1,0,0) (r.l.u.). The zero-field energy is given by 
^(Q) — ^[1 ~ 2(7q]2. Consequently the soft mode indicating transition to (zero-field) AFQ order would occur at 
gQ = 0.5 which is larger than the above value of 0.3 for pure PrOs4Sbi2. Therefore application of a magnetic field 
is necessary to achieve a soft mode uj+{Ql) = at a critical field h^.- The dispersions in Eq. (I3ip are approximations 
where the field dependence of the F^ — F° splitting has been neglected. This is possible as long as their dipolar 



9 



matrix element (P,S^ 1 which is true for the case cP = 0.067 (Fig. [2] (left)) Then the soft mode condition leads to 
the approximate critical field he/ A = (1 — 2gQ)2 above which AFQ order will be induced. Using gq = 0.3 leads to 
he/ A = 0.586 which is close to the exact value 0.632 given in Ref. [l^. 

Calculation of electron self energy and mass enhancement in the dispersive case proceeds now exactly along the 
lines described in Sect. IIIII The main modifications arise from the fact that a more sophisticated approximation for 
the electron-quadrupolar cxciton spectral function is employed. For large cut-off energies ec oo one obtains 




which closely parallels the expression derived by Nakajima and Watabe for the effective mass enhancement due 
to electron-phonon interaction. Here we use q = qq with q = 2pF cos 6 (Fig.[T]) where q has polar and azimuthal angles 
9 and 4>, respectively. Due to the geometric restrictions only half the solid angle {2tt) contributes in the momentum 
integral. Replacing the dispersive modes ti;„(q) by the dispersionless singlet-triplet excitation energies (5„ leads to the 
previous result in Eqs. (|21l25p . Using the explicit matrix elements and dispersions we may again represent the mass 
enhancement in the form of Eq. (PT|) 

^ = 1 + 9'JV(0)/w(H) 

Here the first and second terms are due to the virtual excitations of lo± (q) and loq (q) bosons respectively. In the 
dispersionless limit this reduces to the previous result in Eq. ((25|) for the case (P — (no tetragonal CEF) which 
corresponds to the present treatment due to the neglect of the Eg , level repulsion implied in the dispersions 
of Eq. The above expression for the mass enhancement have to be evaluated numerically due to the BZ 

integrations. This will be discussed in the next section. But we may nevertheless gain some qualitative insights by 
simple approximations to these integrals in the zero-field case. For that purpose we expand aj(q) around one of the six 
equivalent zone boundary X- points with Q =(±1,0,0) etc.. Then one obtains an isotropic approximate dispersion 
given by 

Lu{cir=^Q+^lM' (34) 

where q' = Q — q is the momentum vector counted from the X-point and q' is its length. Furthermore wq — 
A[l — 2(7q]^ = he and luq — y/ggA. On approaching the critical ffg = 5 the soft mode frequency ujq vanishes. In this 
limit the integral in Eq. ([5^ may easily be evaluated. The approximate Fermi surface geometry of PrOs4Sbi2 shows 
that pf — ^ (r.l.u.) or 2pp = \f2 (Fig. [T]). Therefore 2pp > Q which means that the minimum in a;(q) is included 

in the domain of the momentum integrals in Eqs. (I32l33p . We restrict the latter to the sphere around the minimum 
at Q (X-point) with a cutoff radius given by < 1 where only 1/2 contributes due to geometric restrictions. We 
then obtain from Eq. [31] 

^^l + ff^iV(0)(^)^ (35) 
TO 47r cjQ 

where the momentum cutoff around X is defined such that the quadratic expansion in Eq. p4p is still valid. 
Note that although the exciton energy becomes soft ljq — > at the zone boundary there is no divergence in the 
mass renormalisation. This is due to the small phase space volume around the X-point which gives only a small 
contribution despite the vanishing exciton frequency. In addition the singular contribution is suppressed by the fact 
that the form factor 9^9^ vanishes exactly at the X-point directions and only contributions from its environment are 
picked up by the integration. The above expression is formally quite similar to the dispersionless result of Eq. (|27p . 
In the latter case the renormalisation diverges when A — > because this corresponds to a softening in the whole BZ. 
Thus we conclude from Eqs. (|27l35p that the inclusion of a mode dispersion removes the problem of singular mass 
renormalisation, m* /m stays finite for all coupling constants gq, even when gq = gq when the exciton frequency 
becomes soft at Q. However the above formula cannot give a reliable estimate for to* /to due to the strong momentum 
cutoff dependence, we therefore have to employ a numerical evaluation of Eq. (|33p . 
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FIG. 5: Quadruplar susceptibility Axq ~ 5m* /m as function of effective intersite coupling gg. For zero field the intersite 
coupling or dispersive width of uj{ci) has little effect on the mass enhancement. It actually decreases slightly when Qq — 0.5 is 
approached despite the appearance of the soft mode tj(Q) — ^ 0. At the critical field he the dispersion has much stronger influence: 
when the latter becomes small (decreasing qq) the mass enhancement at he strongly increases. For PrOs4Sbi2 {gg — 0.3) one 
may expect little field dependence of Sm* /m between h = and h = he- The inset shows the dependence of the AFQ critical 
field he on gg with the approximation ~ 0. 

VI. NUMERICAL RESULTS AND DISCUSSION 

We first discuss the mass enhancement in the large band width approximation ec ^ cxd for dispcrsionlcss undamped 
CEF excitations. The absolute value of m*{h)/m = 1 + 6m* /m is determined by the effective coupling constant 
g = g'^N{0). Approximating A^(Q) 1/2W (2W=band width) this may be written as g/A = ^^{W/A). Here we 
introduced the dimensionless quadrupolar coupling constant A — gN{0). Assuming typical values of W =1 eV, 
A ~ 0.02 and using A = 8 K we obtain g/Af ~ 0.077 as the size of the effective coupling for the quadrupolar mass 
enhancement mechanism. Using cP = 0.067 and hence ~ 32.6 we obtain a zero field enhancement of m* /m ~ 16. 
This corresponds to the right magnitude for thermal mass enhancement. 

For discussion of the field dependence we use the quadrupolar susceptibility which contains only d as adjustable 
parameter to avoid specifying g. In the case of weak tetrahedral CEF such as realised in Pri_2:Laa;0s4Sbi2 the 
level repulsion of Fg and r° is also weak and therefore the F^ level crosses the Ts ground state at a critical field 
he = A/(l — 5^) in the dispersionless case. The decrease in the excitation gap for h < he and the field dependence of 
matrix elements leads to a field dependence of XqW (Eq- (EH)) which is shown in Fig. [2] for = and (P = 0.067. 
For larger cP the increase is diminished and eventualUy for (P > 0.42 the level repulsion due to the tetrahedral CEF 
is strong enough to lead to an increase in excitation energy and hence to a decreasing effective mass. Likewise the 
level repulsion prevents field induced AFQ order, therefore it is not appropriate for concentrated PrOs4Sbi2. This 
case resembles more that of pure metallic Pr where a singlet-singlet level repulsion in a field also leads to a decrease 
in quasiparticle mass [l3| . 

Let us next turn to the divergent mass renormalization which is predicted for dispersionless undamped CEF excita- 
tions when the triplet level E:^ approaches the singlet ground state level Eg (Fig.[2|). Therefore in the (level crossing) 
case < 0.42 this approach is only valid for moderate fields. As the divergence follows directly from the general 
analytic structure of the corresponding electron self energy it persists also in the self consistent solution for finite 
band width. Selfconsistency leads to an overall reduction of the mass renormalization as can be seen from Fig. 3) 

The (unphysical) divergence of the mass enhancement close to the critical magnetic field is an artefact of the model 
which assumes dispersionless undamped CEF excitations. Inelastic neutron scattering [2l[ however, have shown that 
the singlet-triplet excitations have a pronounced dispersion. The dispersive width corresponds to ~ 40% of the CEF 
excitation energy A. An applied field of critical strength therefore leads to a softening of u'(q) only in the restricted 
phase space around the AFQ ordering vector Q. Consequently the mass renormalisation will be finite even at the 
critical field he for AFQ order when ci;(Q)=0. This is shown in Fig. [3^ for various effective quadrupolar coupling 
strenghts gg. For gg = the mass renormalisation at he would diverge as in Fig. [5] because the excitation energy 
becomes soft for all q-vectors. For small gg and dispersive width the softening apperars only around Q (Fig. [8)3) 
and the mass enhancement is finite, though still large at he- It is progressively diminished with further increasing gg 
because the softening becomes strongly constricted around Q (see Fig. [SJj). For the value gq = 0.3 appropriate for 
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PrOs4Sbi2 the dispersion is sufRciently large to suppress the field dependence of 5m* /m as depicted in Fig. [3^ (full 
line). 

An alternative presentation of results for the dispersive case is shown in Fig. [5] The mass enhancement is shown 
as function of quadrupolar coupling gg (~ dispersive width) for the two limiting cases h — and h = hc{gQ). For 
large gq and dispersion the field variation of dm* /m between h ^ and h = he becomes small. This is partly due 
to the fact that hcigg) itself becomes small for large gg (see inset of Fig. [SJ. When gg decreases the difference in 
dm* /m for h = Q,hc increases rapidly because the mass enhancement at h = he becomes singular when approaching 
the dispersionless case gg 0. The arrow corresponds to the proper value of gg for PrOs4Sbi2 and it shows again 
that one should expect little field dependence of the mass enhancement in this case. 

VII. CONCLUSION AND OUTLOOK 

In this work we have studied in detail the quasiparticle mass enhancement originating in the aspherical Coulomb 
scattering of conduction electrons from singlet triplet CEF excitations. This model has some relevance for the heavy 
fermion superconductor PrOs4Sbi2 where the Pr'^+ 4f states are subject to a tetrahedral CEF leading to a singlet 
ground state and an excited triplet. For small tetrahedral CEF the latter has a mostly nonmagnetic character and 
therefore may be excited by aspherical Coulomb scattering from conduction electrons. These virtual second order 
processes lead to a quasiparticle mass renormalisation which may well be the source of the large thermal and dHvA 
effective masses observed in PrOs4Sbi2. A hybridsation mechanism between conduction and 4f electrons can be ruled 
out since the Fermi surface of PrOs4Sbi2 is identical to that of LaOs4Sbi2 which advocates for fully locallised 4f 
electrons in Pr. Indeed the well defined CEF excitations seen in INS [2l| support this view. 

If aspherical Coulomb scattering of conduction electrons plays a role in the mass enhancement one should expect a 
field dependence of the latter because the triplet excited state splits in the field. For small enough tetrahedral CEF 
characterised by the parameter <C 1 the lowest triplet component crosses the singlet ground state at a critical field 
he- The mass enhancement in second order perturbation theory with respect to Hac then increases with field and 
becomes singular at he- For larger tetrahedral CEF (d^ > 0.42) the excitation energy between singlet ground state 
increases with field leading to a decrease of the mass enhancement, similar as has been observed in Pr metal where 
the mass renormalisation is due to exchange scattering from a singlet-singlet CEF level scheme. 

The observed singular mass enhancement close to the critical field of level crossing is an artefact of the dispersionless 
model, both in the perturbative and selfconsistent treatment. Any dispersion of the singlet-triplet excitations due to 
effective intersite quadrupolar interactions will lead to a finite effective quasiparticle mass. We have shown that the 
enhancement decreases strongly with increasing dispersion because the phase space for conduction electron scattering 
from low lying CEF excitations (quadrupolar excitons) is constrained to the wave vector Q of incipient field induced 
AFQ order. For a quadrupolar coupling constant gg — 0.3 corresponding to PrOs4Sbi2 the field dependence is 
reduced to a few percent. In addition this compound is superconducting below Tc = 1.85 K with Hc2 = 2.2 T and 
has a field induced AFQ phase above He = 4.5 T. Therefore only a reduced field range is left to observe the small 
field dependence possible at gg — 0.3. We conclude that the concentrated PrOs4Sbi2 is not a favorable system to 
observe the field dependent mass enhancement due to aspherical Coulomb scattering. 

A more promising system may be the La-diluted systems Pri_a;La2:Os4Sbi2. On increasing x the average distance 
between the Pr 4f shells becomes larger and therefore the effective quadruplar coupling gg{x) will decrease with 
X, i.e. the dispersion of the 4f CEF ecxitons will progressively decrease. This means that the field dependence of 
effective masses will become more pronounced according to Figs. I3I5I Of course the absolute (zero-field) size of the 
mass enhancement is also reduced since the self energy in Eq. (1171) will be proportional to the number (1-x) of Pr 
sites. There should however be an intermediate concentration region for x where the field dependence is pronounced 
(gg small) and the ^{x) still large enough as compared to the other (lattice or CEF-Schottky) contributions such 
that the field dependence of 7(2;, H) is experimentally accessible. Furthermore in this region of x one may probe a 
larger field range because there is no more AFQ order present. Therefore we propose that the field dependence of the 
electronic specific heat in mixed crystals of Pri_a;La2:0s4Sbi2 is systematically investigated and analysed. It may 
hold important clues to the microscopic origin of the large effective mass in the concentrated compound PrOs4Sbi2. 
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